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\^ Abstract. We have carried out a sample survey in UBVR 

^\ photometry and proper motions in different directions in 
^ the Galaxy, as part of an investigation of galactic struc- 
^ ture and evolution. Three fields in the direction of galactic 
anticentre (/ = 167°, b = 47°), galactic centre (/ = 3°, 
b = 47°) and galactic antirotation (I = 278°, b = 47°) 
have been surveyed. Using our photographic photometry, 
we determine photometric distances for a sub-sample of 
stars in the color range 0.3<(B-V)<0.9. The stellar space 

^ velocities (U, V and W) are derived directly from the mea- 
sured proper motions and distances. 

^-H Using our new data together with wide-area surveys in 
other fields available to date, we discuss the radial and ver- 



tical structure of the Galaxy. We have derived the density 



in 

laws for stars as a function of distance from the galactic 
i-G plane for each absolute magnitude interval. The density 
^^laws for stars with My >3.5 follow a sum of two expo- 
^ nentials with scale heights of 260±50 pc (thin disk) and 
^ 760±50 pc, respectively. This second exponential corre- 
^ sponds to a thick disk with a local density of 7A't_\'l % 



relative to the thin disk. The scale lengths for these two 
populations are respectively 2.3±0.6 kpc and 3±1 kpc. 

The kinematical distribution of F and G-type stars 
have been probed to distances up to 3.5 kpc above the 
galactic plane. A new value for the solar motion has been 
determined from moderately distant stars (l<z<2 kpc). 
It is consistent with local determinations and implies that 
there is no large motion of the LSR relative to the mean 
motion of stars at 1-2 kpc above the galactic plane. The 
rotational velocity curve is found flat in the solar neigh- 
borhood. The radial gradient in velocity dispersions has 
been determined for the thin disk population. The thick 
disk appears as a kinematically distinct population from 
the thin disk and shows no vertical gradient. A multivari- 
ate discriminant analysis is also used to distingiush the 
thick disk from the thin disk and to estimate its asym- 



metric drift. It is found to be 53±10 km/s, independent 
of the galactic radius. Of the many models that have been 
proposed for the origin of the thick disk, the evidence at 
present seems to favour a model in which thick disk formed 
through the rapid dynamical heating of an early disk by 
satellite accretion into the disk. 
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Stellar content - Galaxy: structure - Methods: statistical 
- Surveys 
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1. Introduction 

The concept of stellar populations has played a critical 
role in the development of modern studies of stellar evo- 
lution, galactic structure and galactic evolution. To un- 
derstand fully the stellar populations and their differences 
would take us a long way toward a complete theory of 
galactic structure and evolution. Yet, there are many dif- 
ficult questions that have been very unyielding to final 
solutions, such as, how can we define the populations in 
a way that is physically meaningful, how many discrete 
populations really are there, and how can we unambigu- 
ously separate the stellar constituents of the populations 
for detailed analyses ? 

It was only with the discovery of the asymmetric drift 
of high- velocity stars (Stromberg 1924; Oort 1922, 1926), 
it became clear that multiple overlapping components 
would be needed if an ellipsoidal model were to provide an 
adequate representation of the data. After the identifica- 
tion of the galactic centre (Shapley 1918), and the discov- 
ery of the differential rotation of the galactic disk (Oort 
1922, 1927; Lindblad 1925, 1927, 1959; Strfimberg 1924), 
the multicomponent concept was taken to the limit with 
Lindblad's 1927 (see also 1959) proposed division of the 
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Galaxy "into an infinite number of 'sub-systems' of vary- 
ing mean speed of rotation". Stromborg's (1924) rompro- 
hensive analysis of the available data on stellar distances, 
proper motions, and radial velocities in the mid-1920s was 
among the first to show the differences in the velocity el- 
lipsoids of various stellar groups. His Icinematic separation 
into groups was the beginning of the classification of the 
galactic components that eventually led to the population 
concept. 

The pioneering work on populations was carried out 
by Baade during the 1940's and 1950's using mostly ob- 
servations of galaxies in the Local Group. In 1944, Baade 
introduced the idea of two basic stellar populations (Baade 
1944, 1958a). Unlike the kinematic components discussed 
previously by Stromberg (1924) and Oort (1922, 1926), 
Baade's two populations were defined according to the 
morphology of the rospoctivo color-magnitude diagrams. 
Population I stars are stars with pliotometric character- 
istics similar to those found in open clusters and stel- 
lar associations (e.g. luminous early-type stars, classical 
Cepheids), and Population II stars are similar to the stars 
found in globular clusters (e.g. RR Lyrae stars, subdwarfs. 
Population II Cepheids). 

At the conclusion of the 1957 Vatican conference on 
stellar populations, five stellar populations were advocated 
: extreme population I. older population I. disk popu- 
lation, intermediate population II, and halo population 
II. Excellent reviews on the subject have been written 
by King (1971). van den Bergh (1975), Mould (1982), 
Sandage (1986) and others. 

Throughout the last decade there has been further de- 
bate over the structure and history of our Galaxy. One 
of the main controversies concerns the existence of a thick 
disk, i.e., a component with characteristics intermediate to 
those of the thin disk and halo populations, and whether 
or not three components (rather than two, five or more 
components) are needed to account for star counts, proper 
motion distributions and other observations. There are 
now strong evidences that our Galaxy is not simply de- 
scribed by the traditional disk and a halo, but that there 
must exist a third component, a thick disk, which is dom- 
inated over distances on the order of kiloparsecs above 
the plane, and intermediate in kinematics and metallic- 
ity between the disk and halo populations. The idea of 
"third component'' lias apparently been around in one 
form or another for decades, at least since the time of 
the 1957 Vatican conference (see Blaauw 1965). The con- 
cept of a thick disk came into existence after Gilmore & 
Reid (1983) showed need for it through their inability to fit 
two-component models to starcount data. Since Gilmore 
& Reid (1983), other starcount surveys have reached sim- 
ilar conclusions (Yoshii et al. 1987; Sandage 1987; Fenkart 
1988). Proponents of a three-component model have raised 
questions about the nature of the thick disk, in particu- 
lar, whether it is kinematically discrete (Wyse & Gilmore 
1986; Sandage & Fonts 1987; Gilmore, Wyse & Kuijken 



1989; Soubiran 1993ab) and/or chemically discrete (Car- 
ney et al. 1989: Carney. Latham k Laird 1990). Majewski 
(1992) have shown that an intermediate component must 
exist in order to account for a bimodal distribution in 
both the V-velocity and UV-excess up to 5.5 kpc above 
the galactic plane, far beyond the region dominated by 
the traditional thin disk of scale height ~ 260 pc or less. 
The population of "disk" globular clusters and the appar- 
ent excess of intermediate - metallicity stars with vertical 
scale height h^ ~ 1 kpc, represent also the evidence for 
a thick disk. (Baade 1958b; Morgan 1959; Zinn 1985; Ar- 
mandroff & Zinn 1988). 

It should be stressed here that starcount analysis alone 
do not warrant the reality of thick disk stars as a third ma- 
jor component of the Galaxy, because adding new compo- 
nent (s) in models always gives a superior fit to the data. In 
fact, current interpretations of thick disk stars have been 
ranging between a separate component (Gilmore & Wyse 
1985; Bienayme et al. 1990, 1992) and a high-velocity tail 
of the old disk (Norris 1987ab; see also Norris & Green 
1989). This third discrete population could be the signa- 
ture of a merger event early in the Galaxy's history oc- 
curing shortly after the disk had formed (Carney et al. 
1989). In fact, the arguments regarding the existence of 
this third population bear on problems of galaxy forma- 
tion and evolution. For example, "the G dwarf problem", 
that the local number of metal-poor long-lived stars is 
too small to be consistent with simple models of chemical 
evolution, may be solved straightforwardly by use of an 
intermediate-metallicity "reservoir'' dispersed above and 
below the plane (Gilmore & Wyse 1986). 

Studies of stellar populations often take as their start- 
ing point large-scale surveys, with selection criteria involv- 
ing one or more observable variables. The definition of the 
survey has to be defined in order to ensure a good retrieval 
of the population parameters we look for. The physical lo- 
cation of the surveys is arguably the most important con- 
straint. Obviously several survey locations need to be cho- 
sen in order to vary the mixing proportions of the various 
populations and to determine the scale heights, lengths 
and relative densities. On the other hand, the survey must 
probe to great distances so that the regions dominated by 
each of the Galaxy's populations are within the limits of 
the survey. For example it is necessary to reach V magni- 
tude about 19 in order to reach domains dominated solely 
by the halo populations. Finally, in order to be confident 
that the data are free of systematic bias, it must l)e com- 
plete, i.e. all stars within a region of the sky must be 
surveyed without any selection criteria which bear on the 
problems being studied. 

With these points in mind, we have carried out a sam- 
ple survey in UBVR photometry and proper motions in 
3 directions in the Galaxy. The 3 fields chosen are in the 
direction of galactic anticentre (/ = 167°, b = 47°; Ojha et 
al. 1994a, 1995; hereafter GAG 1,2), galactic centre (/ = 3°, 
b = 47°; Ojha et al. 1994b; hereafter GC) and galactic an- 
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tirotation (/ = 278°, b = 47°; Ojha et al. 1995; hereafter 
GAR). In this paper, using our new data (GAC1.2. GC 
& GAR) together with similar data set towards the poles, 
we derive the kinematical and structural parameters of the 
thin disk and thick disk components of the Galaxy. 

In the present study, we shall take the view that within 
a particular volume of the Galaxy, one deals with a fi- 
nite mixture of discrete stellar populations. Our approach 
has been to see how far one can go using gaussian mix- 
tures to model stellar populations. Thus, our discussion 
of mixture models will be limited to the following vari- 
ables : the three kinematical variables (U, V & W), metal 
abundance [Fe/H] or distance. Under the assumption that 
within a stellar population, each of these variables fol- 
low a gaussian distribution, overall mixture distributions, 
and posterior mixing proportions are computed for two 
or three-component models. Parameters and error esti- 
mations are discussed in general terms. This (minimum) 
system of variables provides an ample basis for the dis- 
cussion of issues that are relevant to the study of stellar 
populations. The purpose of the present work is : 

(i) To present a statistical approach to the analy- 
sis of stellar populations through the application of fi- 
nite mixture model (Stochastic-Estimation-Maximization 
method). 

(ii) To examine the statistical properties of multivari- 
ate finite mixture models. 

(iii) To estimate the parameters that describe the un- 
derlying stellar populations and for determining the num- 
ber of discrete stellar populations. 

The outline of the paper is as follows: In §2, we discuss 
briefly the use of multicolor photometry to derive the ab- 
solute magnitude and stellar distance for each star in our 
surveys. The derivation of the stellar space velocity UVW 
is explained in §,3. In §4, we present the variation of the 
stellar number density as a function of distance from the 
galactic plane. From this, we deduce the structural param- 
eters of the thin disk and thick disk populations. In §5. we 
derive the kinematics for the thin and thick disks from 
two statistical methods. In §6, we derive the components 
of the solar motion and discuss the consistency between 
the kinematics and density scales of the two populations 
(thin and thick disks) via the asymmetric drift relation. In 
§7, we discuss briefly our main results and the thick disk 
formation scenarios. 

2. Stellar distance estimates 

Photometric parallaxes were determined by estimating ab- 
solute stellar magnitudes, which were obtained from a My 
versus B-V relation and taking into account the metallic- 
ity change as a function of the distance from the galactic 
plane. We made use of the BV photographic photometry 
to derive distances for F and G-type stars in our samples. 
The absolute magnitude-color relation for main-sequence 
stars of solar metallicity is derived as : 



Mv = 5.83(5- y) 0.7.53 

This relation is valid on -0.1<B-Y<1.5 color index 
range. A similar color-magnitude relation has been cal- 
culated for the Hyades main-sequence stars (Reid 1993). 
While there is good agreement between the two relations, 
the Hyades relation contains higher order of terms in B-V. 
Since the three fields (GAC1,2, GC & GAR) are at inter- 
mediate latitudes (|6|=47°), therefore we can neglect the 
influence of interstellar extinction in distance calculation. 

2.1. Absolute magnitude correction 

The color-absolute magnitude relation changes with 
metallicity so that metal-poor dwarf stars are fainter than 
metal-rich stars of the same temperature or color, and the 
opposite is true for giants. Consequently, due to the ver- 
tical metallicity gradient, the distance or absolute magni- 
tude derived must be corrected according to the distance 
from the galactic plane. 

The local value of the vertical metallicity gradient 
d\Fe/H] J ^ y determined. Current estimates cover 

dz 

most values from -0.225 kpc^^ (.lanes 1979) to -0.8 
kpc~^ (McClure & Crawford 1971; Jennens 1975; Yoss 
& Hartkopf 1979). Recent studies of field stars RR Lyrae 
(Butler et al. 1979, 1982) and globular cluster (Zinn 1980; 
da Costa, Ortolani & Mould 1982) abundances show lit- 
tle if any systematic trend in [Fe/H] within 1-2 kpc of 
the galactic plane, a dominance of ol)jects with [Fe/H] ~ 
-1.5 at 5 kpc, but little or no further systematic decrease 
beyond this distance. 

This evidence for a smooth gradient from 1 to 

5 kpc is complicated by the results of Rodgers, Harding 

6 Sadler (1981) who derived solar abundances for main- 
sequence A stars up to 4.5 kpc from the plane, and of 
Yoss and Hartkopf (1979) who also derived near solar CN 
abundances for K giants up to 3 kpc. 

In view of these uncertainties, Kuijken & Gilmore 
(1989) produced a model of the vertical metallicity gra- 
dient as indicative rather than conclusive estimates of 
the true gradient. Following Kuijken & Gilmore (1989) 
or Soubiran (1992), we assume the vertical gradient of 
metallicity is : 

^^^^^^"^ = _0.3 kpc-^ z < 5400 pc 
dz 

[Fe/H] = -1.5 z > 5400 pc 

with [Fe/H] = +0.12 for z = 0. 

The most recent study of the sensitivity of My to 
[Fe/H] is by Laird et al. (1988). From a study of local sub- 
dwarfs with trigonometric parallaxes, they find the follow- 
ing relations between the metallicity-dependent UV-excess 
photometric parameter ^o.g and the absolute magnitude 
difference AMy, relative to a Hyades main-sequence star 
of the same B-V color : 
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AMv = ( 2.31 1.04(i? ^) ^(_o.6888^o.a + 53.1<a 
l.oy4 

-97.004^jj g) < ^0.8 < 0.25 

and 

AMv = ,5o.6 < 

The UV excess - metallicity relation was studied by 
Carney (1979), and can be approximated as : 

^o g = -0.0776+ VO.OllOl- 0.05353[Ffi/iJ] 

So we obtain a AMy correction that depends on z. 
The absolute magnitudes calculated from the My versus 
B-V relation was thoroforo corroctod by these amounts. 
In practice, a new distance is calculated, a revised AMy 
correction applied, and the procedure was repeated to con- 
vergence. 

2.2. Errors in lUstance w.easurements 

A variety of systematic errors affect the determination of 
stellar distances. The first source of errors could be from 
B and V magnitudes. This affects the distances of individ- 
ual stars, but should not affect mean distance of a sample 
of stars, since the photometric errors should be randomly 
distributed. More important are the systematic errors that 
could be produced by the contamination by subgiant, gi- 
ant and binary stars in our samples. 

2.2.1. Subgiant contamination 

We used the main-sequence color-absolute magnitude re- 
lation to derive the distances. For the distance determi- 
nation, it is firstly assumed that all stars are unevolved. 
This is certainly not applicable to all stars in our survey, 
but without spectroscopy or photometric gravity indica- 
tors there is no a priori way to separate main-sequence 
dwarfs from evolved subgiant stars at the main-sequence 
turnoff. In our selection of subsamples of F and G-type 
stars (3.5<My<6 or 0.3<B-A'<0.9). the contamination by 
giant stars should be small because at the above absolute 
magnitude and color intervals, distances and velocities im- 
plied for red giant stars are improbable. However contam- 
ination of main sequence (4<My<6) by subgiants with 
2<My<4 (with same color) could cause a serious under- 
estimate of distance. 

We have made a selection by color, by raising the 
blue limit of the color range (B-V<0.9) beyond the re- 
gion where subgiant stars are expected to contribute sub- 
stantially and we can greatly reduce their contamination. 
Only drawback of this selection is that we loss the farthest 
distance and most of the halo stars. The Besangon model 
of population synthesis (Robin & Creze 1986; Bienayme 
et al. 1987) has been used to check the proportion of sub- 
giants in our subsample of stars. If we suppose that there 



exists an intermediate-metallicity population of stars with 
scale height li;~750 pc and local density of 6% relative to 
the thin disk, then we would expect 30% of the intermedi- 
ate subgiant population (My~4) to have (14<V<16) and 
the bulk to have V'~17. An old population I distribution 
(h^~250 pc) would have subgiants (My'~4) contributing 
mostly at V<14.5. The number of possible disk and thick 
disk subgiants in our survey is minimal, since the total 
number of stars in the appropriate magnitude and color 
range is about a 10%. 

For the halo, the distribution is fairly broad. We would 
expect the largest contribution from halo subgiants at 
15.5<V<19.5. The worst contamination of the main se- 
quence by spheroid giants probably occurs at 0.6<B- 
V<0.8 or 5.9<My<6.7 (Sandage 1982, Bahcall et al. 
1983). Recently, a spectroscopic observations have been 
performed on a sample of stars in a field of 9.6 scjuare 
degrees in the galactic centre direction (/ c± 3°, h 45°) 
by Perrin et al. (1995). The most striking result of these 
observations is the very low proportion of Population II 
giants in this direction, about 2 times less than the pre- 
dictions of the current galactic models. Similar result has 
been found in other intermediate latitude fields by differ- 
ent authors (Friel 1987, Morrison 1993). 

It is possible that some ranges of (V, My) are not 
seriously contaminated. For example, if in the range 
4<My<6, the luminosities of the subgiants are close to 
their main sequence values then the mis-identified evolved 
stars do not come from a much larger volume (Tritton & 
Morton 1984). Also at distances where the disk density 
dominates over the spheroid, subgiant or giant contami- 
nation should not be a problem. 

2.2.2. Binary contamination 

Another systematic error affecting the distance determi- 
nation is the problem of unresolved binaries. Every star 
has been treated as if it were a single object, but this 
is certainly not the case for a significant portion of the 
sample "but unknown". Only binaries with components 
of approximately equal mass are problematical. In this 
case, the distance error will be 40%, while for slightly un- 
equal mass stars, there is a concern for color errors in 
the primary (a 10% error in color results in 0.6 mag er- 
ror in estimated absolute magnitude). The binary fraction 
of Population I stars is supposed to be relatively large- 
approximately 55%. The thick disk binary fraction is es- 
sentially unknown, but it is likely to be somewhere be- 
tween that of the thin disk and the halo. In summary, sys- 
tematic errors in the distance determination due to the 
presence of binaries are probably most severe for young 
disk stars. Clearly this problem would be solved with a 
systematic spectroscopic survey for the sample of stars at 
moderately high resolution to identify spectroscopic bina- 
ries. Recently, such a programme has been initiated by 
Soubiran (1994) for subsamples of stars in our surveys. 
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2.2.3. Photometric errors 

The major factors contributing to distance errors are the 
photometric errors in the color ctb-v au'l magnitude ay. 
To roughly estimate the contribution from all of these, it 
is necessary to go back to the observable quantities, as 
some of these errors are correlated, e.g. rrD_y and (Tjpf/f;] 
with omv ■ From our in situ surveys, we expect the mean 
errors in (Jb-v = 0.1 and ay = 0.07 for V = 11 to 18. 
From these values, our estimate of error to be about 20 % 
in distance seems to be realistic up to z = 3 kpc. However, 
we have estimated that for farthest distances (z > 3 kpc), 
there may be an error in our distances as large as 30 %. 

3. Tangential velocities 

The cardinal components of the stellar space velocity (in 
km/s). U, V and W (where U is defined as positive in the 
direction of the Galactic anticentre. V is positive in the 
direction of Galactic rotation, and W is positive in the 
direction of NGP) were derived from proper motions, /.// 
and /ij (in arcsec year^^). and distance d (in pc) for the 
three fields as follows (Murray 1983) : 

U\ / 4.74 dfii' 

V ] =Rg\ 4.74 d in 

Wj \ Vr 

The matrix Rg for each field is as follows 



/ 0.7497 0.0368 -0.6608\ 
GACl,2 0.0808 0.9859 0.1466 

V 0.6569 -0.1633 0.7381 / 

/ 0.5301 -0.5049 0.6813\ 
GC 0.6674 0.7440 0.0320 
\ -0.5230 0.4377 0.7314/ 

/ 0.8753 -0.4746 0.0931 \ 
GAR 0.4043 0.6123 -0.6794 
\ 0.2655 0.6323 0.7278 / 

The unknown radial velocity contributions (Vr) to 
UVW are as follows and are neglected in the following 
study : 



GAC1,2 : Uy^ 
0.738 Yr 



0.661 Vr; V 



Vr 



0.147 Vr; W 



Vr 



GC : = -0.681 Vr; Vy, = 0.032 Vr; Wy, = 

0.731 Vr 

GAR : Uy^ = -0.093 Vr; Vy^ = -0.680 Vr; Wy^ = 
0.728 Vr 

It is clear from the above 3 matrices that we are mea- 
suring the velocities (U-W,V) in the direction of GAC1,2, 



(U+W,V) in the direction of GC and (U,V-W) in the di- 
rection of GAR field. We should point out here that the 
measured proper motions are directly converted to U ±W 
and V velocities. For GC or GAC1,2 field : 



U ±W 
and we define : 



4.74 d fib 



id y ~ 4.74 d fii 



4. Space density 

We are interested in deriving the structural parameters 
(e.g. scale height and scale length) of the thin and thick 
disk populations using our data sets at intermediate lati- 
tude. For this we have calculated the stellar space density 
as a function of both absolute magnitude and distance 
from the Galactic plane. The calculation is straightfor- 
ward : 

The space density of total stars (Population I -|- In- 
termediate Population -|- Population II) of the luminosity 
group (Ml < M < M2) falling into the partial volume 
AVi 2 is given by : 



p(f) = p(ri,r2) = 7Vi,2/Ayi 



.1,2 

TVi 2 being the total number of stars of the lumi- 
nosity group in question, and partial volume AVi 2 = 
(7r/180)^(n/3)(r| — rf). ri and r2 are the limiting dis- 
tances; □ : field size in square degrees; f = ((rf -|-r|)/2)^/^ 
: centroi'd distance of AVi 2- 

4.1. Density laws 

The space density distribution can be approximated by a 
double exponential : 



p{R, z) = p{Ro)exp 



R — Rq 



exp 



^ (1) 



Ro = 8.09±0.3 kpc (Pont et al. 1994) is the solar dis- 
tance from the galactic centre. /3(Ro) is the solar "neigh- 
borhood" normalisation. R is the galactocentric distance 
projected upon the galactic plane and z is the height above 
the galactic plane, and Hr are the scale height and scale 
length, respectively. 

The density can be expressed as a function of the dis- 
tance d along the line of sight for a field of galactic coor- 
dinates (/,6) with, 



R= (RI+ d'^cos'^b - 2Rod cosh cost) 



1/2 
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We derived the logarithmic space density functions of 
the thin and thick disk stars in 3.5<My <5 and 5<My <6 
absolute magnitude intervals. In deriving space density, we 
excluded data for all stars with distances outside the in- 
terval 150<z<2750 pc. These precautions are necessary to 
reduce the contamination of FGK dwarf sample by possi- 
ble Population I and II giants, and to keep completeness 
of counts in all studied fields. The 150 pc limit then arises 
both because we have very few stars nearer than V<11 
and because of possible contamination by young disk stars. 
The outer distance limit of z~2750 pc arises because the 
exponential fit is no longer valid beyond this distance due 
to completeness. 

The observed density is fitted by the sum of two ex- 
ponentials of the form shown in Equ. (1). A least-squares 
method is adopted to derive the best model in which the 
exponential distributions of the numl)er density reproduce 
the observed starcounts. Bins of distance have been se- 
lected in order to have at least 50 stars per bin to re- 
duce the Poisson error. The values of the scale lengths of 
thin and thick disk populations determined from the star 
count ratio (see §5.3) were fixed in the calculation. We 
find that varying these scale lengths, the value of scale 
height changes only slightly. To avoid these uncertainties 
in the determination of structural parameters, we have 
used the combination of three intermediate latitude fields 
(GAC1,2, GC & GAR) explained in the next subsection. 
We computed a number of exponential models iteratively 
until the quantity X^^^i [Nf'' - N^'"^^^f is minimized on 
a domain of the major parameters, p is the number of bins. 
It is evident that, in the least-squares solution, the bins 
are weighted nearly proportional to the numl)er of stars 
located in them. We found that the least-squares solution 
for the intermediate component is not greatly changed. In 
figures 1 and 2, we present the density distribution of stars 
for two different absolute magnitude ranges as a function 
of distance above the galactic plane. 

Tables 1 and 2 present the best fitted structural pa- 
rameters for stars with 3.5<My<5 and 5<My<6 abso- 
lute magnitude intervals from the two data sets (GAC1,2 
and GC). 

By comparing table 1 and table 2, we find some dif- 
ferences on scales depending both on fields and absolute 
magnitude intervals. These differences may be due to the 
larger statistical errors in the GAG 1,2 data set (area 7 sq. 
degs. and 5956 stars) compared to the GC data set (area 
15 sq. degs. and 19696 stars). They may also be produced 
by errors in B-V color that is used to select the samples. 



Fig. 1. The density distribution for stars with 3.5<Mv<5 as 
a function of distance above the galactic plane. The continuum 
line represents the sum of two exponentials with scale heights 
given in table 1 and corresponding to the thin disk and thick 
disk 



^.2. The density distnhuiion p(z) for stars with 
fl.5<Mv<5 derived from the three intermediate lat- 
itude fields (GAC1.2. GC and GAR) 

We have used the combination of 3 intermediate latitude 
fields (GAC1,2, GC and GAR) to derive the structural 
parameters of the two populations. To minimize the un- 



Fig. 2. The density distribntion for stars with 5<Mx/<6 as a 
function of distance above the galactic plane. The continuum 
line represents the sum of two exponentials with scale heights 
given in table 2 and corresponding to the thin disk and thick 
disk 
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Table 1. The best fitted structural parameters of the thin and thick disk stars with 3.5<My<5 derived from the two data sets 
(GAC1,2 and GC) 



Field 


Thin Disk 


Thick Disk 


Thin Disk : Thick Disk 


3.5<My<5 


(pc) 


(pc) 


density ratio 


GAC1,2 


336 


807 


100 : 7.3 


GC 


226 


674 


100 : 9.8 



Table 2. The best fitted structural parameters of the thin and thick disk stars with 5<Mv<6 derived from the two data sets 
(GAC1,2 and GC) 



Field 


Thin Disk 


Thick Disk 


Thin Disk : Thick Disk 


5<My<6 


(pc) 


(pc) 


density ratio 


GAC1,2 


267 


849 


100 : 8.6 


GC 


217 


726 


100 : 9.9 



certainty in the determination of density laws, we built 
a more consistent model. We fit simultaneously data in 
each field with the same values of parameters (expect for 
the thin disk scale length that has not been fitted). For 
each disk, we have used the same local normalization of 
density independently of directions. In this case we may 
solve both parameters {Iir & hj) together. The observed 
density is fitted by the sum of two exponentials shown in 
figure 3. The same procedure (described in the previous 
section) was applied to derive the best fitted structural 
parameters. The results are shown in table 3. 

4-3. Discussion 

The best value for the scale height of thin disk is h^ 
= 260±50 pc. In contrast with the generally adopted 
value around 32.5 pc. we find a lower value for the scale 
height. This low value is in agreement with the scale height 
of 249 pc fitted by Kuijken & Gilmore (1989) from a K 
dwarf photometric parallax study in the direction of the 
South Galactic Pole. Haywood (1994, 1995) showed that 
the overall vertical density profile of the galactic disk is 
closed to an exponential with scale height h^~2.50 pc. Ng 
et al. (1995) recent determination gives h^ = 250 pc for 
the scale height of thin disk based on a sample of stars 
toward the galactic centre. 

The thick disk characteristics are, h^ = 760±50 
pc and local density = 7.4lj g % relative to the 
thin disk. Robin et al. (1995) recent determination gives 
h^=760±50 pc, with a local density of 5.6±1.0 % rela- 
tive to the thin disk. This last determination has been 
done with available B, V photometric data taken quite 
extensively from literature. They use a different approach 
to analyse data using a synthetic model reproducing ob- 
servable quantities (magnitude and color counts). This 



Fig. 3. The density distribution for stars (3.5<Mv<5) in three 
fields (GAC1,2, GC and GAR) as a function of distance above 
the galactic plane. The continuum line represents the sum of 
two exponentials with scale height 2C0 pc and 7C0 pc and cor- 
responding to the thin disk and thick disk, respectively 

method is expected to avoid systematic bias that can be 
encountered in inversing the process. Two points can ex- 
plain the small differences between Robin et al. (1995) 
and present determinations. First, the density law used 
by Robin et al. (1995) is slightly different from our own. 
They use a modified exponential with a null derivative at z 
= 0, in order to conform to the potential. Our density law 
is strictly exponential thus the local density deduced from 
the counts at z>l kpc is slightly overestimated in our case 
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Table 3. The best fitted structural parameters for the thin disk and thick disk stars with 3.5<My<5 derived from the three 
data sets (GAC1,2, GC and GAR). The scale length of the thin disk was not well determined due to poor statistics of the data 
in the nearer distance bins (z<800 pc) 



Thin Disk 


Thick Disk 


Thin Disk : Thick Disk 


(pc) 


(pc) hu (pc) 


density ratio 


260±50 


760±50 3800±500 


100 : T-'it'ii 



compared to Robin et al. (1995) result. Second, there is a 
slight correlation between scale height and density when 
determined using star counts (Reid & Majewski, 1993). 
However, Robin et al. (1995) used a larger magnitude in- 
tervals (from V = 12 to 22), thus raising the indetermina- 
tion. 

Our scale height determinations are mainly based on 
the shape of densities on figures 1-3, while the scale lengths 
depend on ratio of counts between opposite fields. Small 
relative photometric errors between fields will not intro- 
duce large change on the shape of figures, but will af- 
fect relative counts between fields and modify the scale 
length. This effect is larger when we select sample with 
small extent in _B — V (e.g. the values of the scale length 
of thick disk range between 2 and 3 kpc for the second 
sample (5<My<6), and 3.7±0.5 kpc from the first sample 
(3.5<My<5)). A conservative estimate of the errors con- 
sists to consider the maximal fluctuations of our various 
determinations for different samples and different fitting 
processes. We found that the scale length is not accurately 
measured but conversly the exact values have no large in- 
fluence on the other measured quantities (scale height and 
local densities). 

We dediice that the scale length of the thick disk 
is liR = 3±1 kpc. Robin et al. (1995) determination 
gives 2.8±0.8 kpc, which is based partly on the same data. 
These results can be compared to recent determination 
given by Soubiran (1993ab), based on a sample of stars 
with proper motions towards the NGP : h^^thm disk='280 
pc. h,j,,,rk r/,.a=700 pc. Von Hippel & Bothun (1993) 
gives h^^thm disk ~290 pc, h^^tMck disk=860±90 pc based 
on a faint Stromgren photometric survey. 

5. Kinematics of stellar populations : 
SEM (Stochastic-Estimation-Maximization) al- 
gorithm 

The kinematical properties of each stellar population are 
related to their spatial distributions. Scale height, velocity 
dispersions and asymmetric drift are linked by the Boltz- 
mann equation. The proportion of each population varies 
with the distance above the galactic plane and the selec- 
tion of a stellar sample at a given distance allows to opti- 
mize the proportion of one population. Since kinematical 
data allow to improve this identification of populations, we 



have minimized bias from mutual contamination of each 
population at a given height by performing a kinemati- 
cal separation of populations. For that purpose we have 
assumed that the kinematic of each population is well ap- 
proximated by a maxwellian velocity distribution. In case 
of an isothermal population and separation of vertical and 
radial motions, this is particularly suited, since no kine- 
matical gradient is expected, mean velocity dispersion and 
asymmetric drift remain constant with height above the 
plane. This allows to recognize one population at various 
heights and to measure its density along the line of sight. 

To perform the kinematical separation, we have used 
a maximum likelihood method (SEM algorithm : Celeux 
& Diebolt 1986) in order to deconvolve the multivariate 
gaussian distributions and estimate the corresponding pa- 
rameters. The aim of the SEM algorithm is to resolve the 
finite mixture density estimation problem under the max- 
imum likelihood approach using a probabilistic teacher 
step. Full details can be found in a review paper published 
by Celeux & Diebolt (1986). Through SEM one can obtain 
the number of components of a gaussian mixture (without 
any assumption on this number), its mean values, disper- 
sions and the percentage of each component with respect 
to the whole sample. This method has already been used 
by Gomez et al. (1990) and Soubiran (1993ab) to charac- 
terize the (U,V,W) parameters of the stellar populations. 

The samples of stars in 2 fields (GAC1.2 and GC) have 
been devided in 6 or 7 bins of distance, and in each bin of 
distance a fit has been performed with a SEM algorithm 
to separate the 2-D gaussian distributions to identify the 
three components (tliin disk, thick disk and halo) of the 
Galaxy. As can be seen in Ojha et al. (1994a) (tables 11 
& 12) and Ojha et al. (1994b) (tables 10 k 11), the per- 
centage ratio of 2 populations (thin disk/thick disk) varies 
with the distance above the galactic plane and therefore 
with the line of sight distance. The thin disk population 
has been identified as a mixture of components (probably 
young and old disks) and thus we observe a gradient in 
velocity dispersions as a function of z. The thick disk pop- 
ulation (table 12 in Ojha et al. (1994a) and table 11 in 
Ojha et al. (1994b)) has been identified as a discrete and 
distinct component. We do not find any gradient either 
radial or vertical in the velocity ellipsoid of the thick disk 
population. This proves that the thick disk is an isother- 
mal population. The halo population has been identified 
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in the farthest distance bins (z>2 kpc) with large velocity 
dispersions and in a prograde rotation. Wo present liere 
new results obtained from the comparison of data from 
different fields that had been studied separately. 

5.1. Kinematics of the thin disk population 

For the small distance intervals r<500 pc, the gaussian fit 
shows that the thin disk-like population was itself divided 
into two components. This may be due to the varying 
contribution of young disk and old disk populations in 
these distance intervals. We observe a continuous increase 
of the velocity dispersion with height. It can be explained 
by the fact that the thin disk has formed over about 10 
Gyr, therefore can be considered as a sum of isothermal 
disks. The dynamical evolution of the orbits results in a 
vertical gradient of velocity dispersion, detected in this 
analysis. 

Using our two data sets in the directions of GC and 
GAC1,2, we obtain a measure of the radial gradient of 

velocitj' dispersion ( ,3 ' ) 1 which is found to be - 
0.18±0.03 kpc^^ or (5.7±1 kpc^^) for stars (in color range 
0.3<B-V<0.9) on a base line of ~ 2 kpc (figure 4). We 

obtain also (^^)=-0.05±0.02 kpc^^. Th ese results are 

in agreement with Fux & Martinet (1994) (-^^ = - 
0.161q'q3 kpc^^). The local value of the velocity disper- 
sion gradient {-^§^) may be determined using the sta- 
bility criterion of the galactic disk introduced by Toomre 
(1964). Based on the Toomre's criterion. Mayor (1974) de- 
duced a value of -0.2 kpc^^ for ^^g^" in contrast to the 
value of zero that one would have in the ellipsoidal hypoth- 
esis. This result has been confirmed by Oblak & Mayor 
(1987) from F, G and K type stars in the Gliese cata- 
logue, who found values between -0.14 and -0.24 kpc~^. 

Our value for is also in good agreement with 

the estimation of Vandervoort (1975) and Erickson (1975) 
obtained from the 3rd and 4th order moments of the lo- 
cal stellar velocity distribution (between -0.19 and -0.23 
kpc^^) based also on marginal stability assumption and 
with Lewis and Freeman (1989) kinematical scale length 

of 4.4±0.3 kpc, (^^ = -0.23±0.02 kpc^^) determined 
using a large number of (about 600) distant old disk K 
giants. 

5.2. Kinematics of the thick disk population 
5.2.1. Gaussian deconvolution 

The SEM solutions of the thick disk derived from each field 
are presented in Ojha et al. (1994ab). The mean kinemat- 
ical results of the thick disk obtained from the 2 fields 
(GAC1,2 and GC) are summarized in table 4. 

Galactic radial gradients of the velocity dispersion 
( "^g"^ ) of the thick disk population were never determined 



Fig. 4. Plot of <Ti7_ as a function of d distances along the line 
of sight at color interval 0.3<B-V<0.9 mag for the two data 
sets : galactic anticentre field (I = 167°, b = 47°) and galactic 
centre field (/ = 3°, 6 = 47°) 

Table 4. Kinematics of the thick disk population (expressed in 
km/s) derived from SEM algorithm. Velocities are with respect 
to the Sun 



Anticentre 


(GAC1,2) 


Centre 


(GC) 


U-W 


-5±6 




4±3 


<V> 


-57±4 


<v> 


-49±3 


'^u^ w 


64±4 




66±3 


(TV 


60±3 


(TV 


56±2 



before. In our present study, the thick disk population ap- 
pears to be isothermal, since the kinematical properties 
are nearly constant along the line of sight. For this reason, 
the velocity characteristics for this population are much 
better defined. We determine a small positive, null and 
negative gradients respectively in GC. CAE and GAC1.2 
fields, that may be explained as a radial gradient of the 
kinematic of the thick disk. Mean values allow to deter- 
mine local properties of the thick disk. From comparison 
of proper motion distributions in each field, the kinemat- 
ical gradients can be estimated on a base line of 5 kpc : 

^^g^'' range from -0.02 to -0.10 kpc~^. 

By combining the kinematical results deduced from 4 
fields (GAC1,2, GC, GAR and NGP (Soubiran 1993ab)), 
we have derived the velocity ellipsoid of the thick disk pop- 
ulation. The mean kinematic parameters are summarized 
in table 5 and discussion is given below. 

Soubiran (1993ab) and Kharchenko et al. (1994) have 
determined kinematical properties of the thick disk from 
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proper motion surveys towards the NGP. Their results can 
be compared to onr results since selection criteria are very 
similar. Applying our distance determination (see §2) and 
SEM to their samples, we obtain nearly the same values 
for dispersions and asymmetric drift. Then we obtain a 
mean ajj = 71±5 km/s from measurements towards the 
NGPs and from our field (GAR). We mention that this 
value is consistent with the dynamical constraint of the 
asymmetric drift relation (see §6). crv had been obtained 
in four fields (three different directions : GC, GAC1,2 and 
NGP) and is apparently constant : mean rry = .57±4 km/s. 

It is difficult to extract aw from our determinations 
of (^Ujiv and ajj using the relation cr^p^ = (afj + (T^)/2, 
since errors from both determinations are added and rriy 
is the (small) difference of two large quantities. The rela- 
tive error on aw is then too large. We prefer to apply an 
indirect determination in order to estimate aw- The scale 
height of an isothermal population is related to its vertical 
velocity dispersion. The exact relation can be determined 
if the vertical potential is known. In fact this method has 
been frequently applied to determine the vertical poten- 
tial from samples of disk stars with known vertical density 
distribution and velocity dispersions. Most recent results 
(Bienayme et al. 1987; Kuijken & Gilmore 1989; Flynn & 
Fuchs 1994) have shown that the vertical potential can be 
explained by the visible mass (stellar and IMS) without 
advocating for a dark matter thin disk. Such determina- 
tions are based on disk stars not very distant from the 
galactic plane, where the vertical force changes linearly, 
and where the scale height of these populations is nearly 
proportional to their velocity dispersions. For thick disk 
stars considered in this paper, their distances range from 
1 to 2 kpc, where the vertical force is nearly constant, 
and so the scale height of the population in this range of 
distance is nearly proportional to a'^y. So aw is less de- 
pendent on the exact value of the potential. Taking the 
potential obtained by Bienayme et al. (1987), Kuijken & 
Gilmore (1989) and Flynn & Fuchs (1994), we determine 
a value of aw= 40 km/s. Assuming the potential of Oort 
(1965) with 30 percent of unseen matter in the disk, we 
would obtain aw = 45 km/s. 

The accuracy on measured dispersions depends on the 
number of stars in the sample. The samples towards the 
NGP are relatively small and result in larger errors (error 
bars given by Soubiran (1993ab) reflect directly the num- 
ber of stars in her sample). Errors on Kharchenko et al. 
(1994) data must be similar (3 to 10 km/s). A second effect 
comes from the accuracy of proper motion that has been 
exactly determined for all these surveys, and range from 
1 to 2 mas/year or 10 to 40 km/s for stars at distance of 1 
to 2 kpc. From the observed proper motion distributions 
and with the help of the Besancon galactic model, we de- 
termine the ellipsoid of the velocity distribution corrected 
from errors on the proper motions (table 5). 

Remark : All previous determinations have been done 
assuming that the main axis of the velocity ellipsoid re- 



mains parallel to the plane. The exact situation is not clear 
while it is sometimes claimed that the ellipsoid should 
point towards the plane in a direction beyond the galac- 
tic centre. If the ellipsoid is oriented in such a way, the 
velocity dispersion observed on the GC field would imply 
a smaller ajj in this direction, and a slightly larger in the 
GAG field. This implies a gradient of ajj, with increas- 
ing values of ajj outwards. This apparently contradicts 
the fact that we observe a small negative gradient for the 
V-dispersions (that is not modified by an ellipsoid incli- 
nation), and since U and V-dispersions are kinematically 
linked. 



Table 5. The mean kinematic parameters of thick disk (in 
km/s) derived from 4 fields (GAC1,2, GC, GAR & NGP). aw is 
determined for the most probable vertical potential (see text). 
Viag is with respect to the Sun 





av 


aw 


^Lag 


Thick disk 67±4 


51±3 


40 


-53±10 



It should be noted that since two years the accuracy of 
the kinematical measurements of the thick disk population 
has greatly improved and the results are well in agreement 
with each others (see table 6). It remains a controversy 
about a possible vertical gradient claimed by Majewski 
(1992) but not seen in data from Soubiran (1993ab), Ojha 
et al (1994ab) and others. The determination of the gradi- 
ent is of great consequence on the scenario of formation for 
this population (see §7). The controversy seems to come 
from the way poptilations are distinguished from each oth- 
ers. The SEM algorithm (used by Soubiran 1993ab and 
Ojha et al. 1994ab) allows to well separate the thick disk 
and to evaluate its kinematics within a good accuracy in 
distance bins where it is preponderant. 

5.2.2. Multivariate discriminant analysis 

Because the kinematics and metallicity of thick disk pop- 
ulation is supposed to differ from the thin disk and halo, it 
may be possible to pick out thick disk stars on the basis of 
the kinematical differences between the three populations. 
We therefore try here to find new constraints on thick disk 
population using samples at intermediate latitude which 
include photometry and proper motions. We have used 
multivariate discriminant analysis (MDA) to qualify the 
thick disk using observations in multidimensional space 
(V, B-V, U-B, III k lit). 

A MDA is used to search for the best discriminant axes 
to project the data in multi-dimensional space to seek the 
optimal separation between different populations. The dis- 
tribution of stars along the discriminant axis (combination 
of the observed axis) shows the clearest separation be- 
tween the thin disk, thick disk and halo. We have applied 
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Autlior(s) 




Thick disk 






<Jv 


<JW 


^Lag 


Wyse & Gilmore 1986 


80 


OU 


OU 


100 


Carney & Latham 1986, 1989 


- 






30 


Norris 1986,1987ab 


- 




60 


20 


Sandage & Fouts 1987 


75 




42 


31 


Ratnatunga & Frficman 1989 


77±16 


54 




— 


Morrison ct al. 1990 


55 


40 




35±10 


Majewski 1992* 


35-130 


35-70 




21-120 


Soiibiran 1993ab 


56d=ll 


43±6 




41±16 


Khardionko ot al. 1994 


67d=l 


54±1 




63±1 


Bartasiuto 1994 


64±5 


49±3 


42±3 


39±5 


Beers & Larson 1994 


63=t7 


42±4 


38±4 


20 


Layden 1994 


59±12 


43±10 


42±10 


34±10 



* The author advocates a vertical kinematical gradient of the intermediate population stars 



this method to our data sets (GC and GAC1,2), using sim- 
ulations from the Besancon model of population synthesis 
to separate the thick disk among other populations and 
to investigate the circular velocity of this population. The 
advantage of this method is that we can extract the kine- 
matic parameters of stellar populations without estimat- 
ing the stellar distances i.e. by using a model of population 
synthesis. 

Using model predictions, we select a subsample of stars 
(where the thick disk stars are in majority), chosen as B- 
V<0.8 and 14<V<15.5. We used the model simulations 
to find the best discriminant axes where to project the 
data in order to separate the thick disk population from 
the thin disk and halo. This was done using a multivari- 
ate discriminant analysis under the MIDAS^ environment. 
Model simulations have been made assuming different cir- 
cular velocities of thick disk. The characteristics of each 
tested model are shown in Ojha et al. (1994b) (table 8). 
To avoid too large Poisson noise in the Monto-Carlo simu- 
lations, we computed at least 10 simulations of 100 square 
degrees for each of the models tested in our analysis. The 
first discriminant axis for the circular velocity of 180 km/s 
of thick disk in the direction of GC is given by : 

X = 0.024(S-y)+0.139([7-S)-0.079y-0.310^/-0.069^6 

It should be noted that the discriminant axis varies 
when we change the circular velocity in model simulations. 
The resulting discriminant axis is dominated by the proper 
motion along the rotation (jii) and by the U-B color index 

^ MIDAS is an acronym for Miinich Image Processing 
and Data Analysis System developed by European Southern 
Observatory. 



due to metallicity differences between the thin disk, thick 
disk and halo. 

To quantitatively estimate the adequacy of the models 
with various circular velocities, we applied a test to 
compare the distribution of the sample on the discriminant 
axis with a set of model predicted distributions assuming 
different circular velocities of thick disk. The is given 
by the following formula : 

2 _ (0; — fl; J 

^ a, 

2=1 ' 

Where n is the number of bins and a h h are the num- 
ber of counts in each bin in the model and observed data 
sets, respectively. The resulting distribution may be 
approximated by a normal distribution given by the ex- 
pression (Wilson & Hilferty 1931) : 

n — i 9(n — 1) 2 

The above expression is approximately normally dis- 
tributed around zero mean with unit variance (Kendall & 
Start 1969). 

Table 7 and figure 5 show the values of the probability 
(in sigmas) of each model to come from the same distri- 
bution as the observed sample. The most probable value 
for the lag or asymmetric drift of thick disk comes out 
to be -55±10 km/s (see figure 5), which shows that the 
circular velocity of thick disk is of the order of 180±10 
km/s (where 'V lsr = 229 km/s and Vp, = 6.3 km/s in 
the model). We notice that the model predictions are at 3 
sigmas of the data. This is due to the fact that the statis- 
tics of the errors in the data is not a Poisson statistics. 
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because of possible systematic errors in the photometry 
and astrometry. 



Table 7. ^^est for models with different circular velocities 
for the thick disk. Lag or asymmetric drift is with respect to 
the Sun, where Vlsjj=229 km/s and V0 = 6.3 km/s in the 
Besangon model 



167°, h = 47°) with best model predictions (with three 
populations) is shown in figure 6. 



V„r (km/s) \Lag (km/s) (in sigmas) 



GC GAC1,2 



150 


-85.3 


8.7 


6.8 


165 


-70.3 


6.8 


5.7 


175 


-60.3 


4.1 


3.8 


180 


-55.3 


4.3 


3.4 


185 


-50.3 


4.9 


3.5 


190 


-45.3 


5.5 


4.0 


215 


-20.3 


6.7 


5.1 



Fig. 6. Distribution over the discriminant axis of observed 
stars (solid points, with -iz^/W Poisson statistics) towards 
galactic anticentre direction; best model predictions (solid full 
thick line) and model predicted stars according to their pop- 
ulations (thick line : thin disk, full thick line : thick disk and 
thin line : halo) 



Fig. 5. test for models with different circular velocities for 
the thick disk population for the two fields (GC and GAC1,2). 
X^ is given in number of sigmas 

We obtain a unique value in both directions, showing 
that no radial gradient seems to occur on a base of 3 kpc 
around the Sun. The distribution over the discriminant 
axis of observed stars (solid points) towards GAC1,2 (/ = 



The resulting asymmetric drift of thick disk (V Lag = - 
55±10 km/s) is in good agreement with the determination 
made using the gaussian doconvolution of velocity distri- 
bution (i.e. \Lag = -53±10 km/s, see §5.2.1). The MDA 
approach we use here to separate the thick disk from other 
populations, allows to avoid two biases. First, by using a 
model of population synthesis, which adopts the same se- 
lections in the model and data, and second, by identifying 
the thick disk from other populations on a physical basis, 
by its metallicity and kinematics and using at the same 
time all 5 observables. 

5.3. Scale lengths from kinematically identified popula- 
tions 

In section 4, we have determined the scale length and scale 
height of the populations from computing the density law 
along the line of sight in different absolute magnitude (or 
color) intervals. The kinematical separation of populations 
using SEM algorithm also gives an estimation of the pro- 
portions and densities of each population along the line of 
sight distance. Figures 7 and 8 show the number of F and 
G type stars of the thin and thick disk respectively, as a 
function of line of sight distance (d/o^). 



D.K.Ojha et al. : Structure and kinematical properties of the Galaxy. 



13 



Using the density function of Equ. (1) the star count 
ratio between tlie two fields GC and GAC1.2 (where \T! — 
Rq\ ~ |2;|), assuming that the luminosity function is the 
same, can be written as : 



AGc{m)/A 



GAC1,2 



hn = 



(m) = exp{+2*\R- Ro\/hR) 
exp{V2 dioa/hn) 

a/2 dtos 



log(AGc('m)/A 



GAC1,2 



(m)) 



By comparing the star count ratio between the two 
data sets in each distance bin, we obtain the scale length 
of thin disk is h^; = 2.6±0.6 kpc and for the thick disk 
hij = 3.6±0.5 kpc. It appears from the test with sim- 
ulated data that SEM results were not very stable when 
two populations have the same proportions. This limit the 
accuracy of density determination in some distance bins of 
transitions. However the results obtained on scale length 
is compatible with the results obtained directly from star 
counts in §4 (table 3), comforting the separation of popu- 
lation made from the SEM method. 



Fig. 7. Tlic observed number of thin disk stars as a function 
of d distances obtained from SEM algorithm from the two data 
sets (GAC1,2 and GC) in 0.3<B-V<0.9 color interval 



6. Solar motion and asymmetric drift relation 
6.1. Solar m.otion 

Velocity components of the solar motion are usually de- 
duced directly from the mean motion of nearby stars rela- 
tive to the Sun. The tangential component V0 of the solar 



Fig. 8. The observed number of thick disk stars as a function 
of d distances obtained from SEM algorithm from the two data 
sets (GAC1,2 and GC) in 0.3<B-V<0.9 color interval 



motion in a circular rotating coordinate frame is based on 
the extrapolation of the asymmetric drift relation (<V>, 
<yfj), where <V> is the mean of the V velocity components 
of a stellar sample, a-j, is the variance. Here, we measure 
the solar motion from distant sample of stars (1 to 3 kpc) 
and compare it to the local determinations. 

Since our samples of stars are non local, we deduce the 
asymmetric drift relation from a general form (see Binney 
& Tremaine, 1987) and we obtain the expression : 



-Vlag = -V4R)+Vrot{R) = 



Ral 

2Vr 



dlnp ^ dlnai 



dR 



dR 



4<> 



1 dp<UW> 



'u 



pdz 



or shortly 



-V^R) + V,ot{R) = 



(2) 



(2.1) 



under the following assumptions that the Galaxy is sta- 
tionary and axisymmetrical, and Viag << Vc. Vc(R) is the 
circular velocity curve and Vrot{R) is the mean rotational 
velocity of stars with radial velocity dispersion ajj. 

The mean apparent velocity of a stellar sample relative 
to the Sun is given by : 

<y >„,,= K„,(i?)-y,(i?Q)-yQ (3) 

combining the equations (2 & 3), we obtain : 
dV^R)^ 



<V>o 



dR 



-{R-Rq)-Vq + 



D 
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In figure 9, we have plotted the (<V>, cr^ ^) points 
(with (T^ ^, = (iTf, + <7^y)/2) obtained from the observa- 
tional data sets in two directions (GC and GAC1,2) for F 
and G-type stars (0.3<B-V<0.9). We find that, for (Tu,w 
< 19 km/s, the (<V>, ct^ ^) points differ significantly 
from the straight line (see also Mayor 1974). In two papers, 
Mayor (1970, 1972) have shown that density waves or local 
perturbations disturb essentially the stellar populations 
with low velocity dispersions. Then discarding these sub- 
populations having a low velocity dispersions and within 
the limits of statistical uncertainties, the extrapolation of 
(7 = gives a mean V0 = 6.6±1.6 km/s. Joined with the 
mean values of the U and W components, we obtain the 
solar motion : Uq = -1.0±6.0 km/s, V0 = 6.6±1.6 km/s 
and W© = 6.4±6.0 km/s. 

The derived standard solar motion is typical for F 
and G-type disk dwarfs. These values can be compared 
to those of Mayor (1974) : (-10.3±1.0, 6.3±0.9, 5.9±0.4) 
km/s and Oblak (1983) : (-8.2±1.8, 5.0±0.7, 5.5±0.4) 
km/s, determined from the local sample of stars in a radius 
of 200 pc around the Sun. There is no significant differ- 
ences between the values of solar motion derived from lo- 
cal and distant samples of stars, implying no visible moan 
motion of local stars relatively to distant one on the merid- 
ional plane. 

6.2. Galactic rotation curve 

The value of D depends on quantities that are equal in 
both opposite fields (GC and GAC1,2). The contribution 
of the gradient of the circular velocity curve has an oppo- 
site sign in the two fields, and comparison of them gives 
a direct measure of Wc(R). On figure 9, the relation 
(< V >, (T^) has nearly the same value for both fields, indi- 
cating a null gradient for the velocity curve. We determine 
VV'c(i?) = -1-0.7 ± 1 km/s based on a range of distances 
of ±1.0 kpc and conclude that the rotation curve is flat 
around the solar radius. Our analysis is in agreement with 
various determinations of the outer rotation curve (Fich 
& Tremaine, 1991) based on more distant tracers. 

6.3. Asymw.etnc drtft 

The gradient of the velocity curve being null, the expres- 
sion of the asymmetric drift can be written in a simpler 
form: 

< V >ots= -Vq + 4 

Simplified expressions for D are frequently found and 
are generally based on assumptions without any strong 
theoretical support. A careful analysis can be found in 
Fux & Martinet (1994). Current assumptions concern the 
last term in the bracket of Equ. (2) that is frequently 
replaced by its value for a spherical or a plane parallel 
potential, but not by its more probable value which may 
be intermediate. For 2 = 0, this last term can be written 
(Fux & Martinet 1994) as ^(1 - ^) with < A < 1. 



Exact value of A depends on unknown quantities like the 
flattening of the dark halo. In the following discussion we 
will choose A = 0.7, according to Fux & Martinet (1994). 

The other assumption concerns the kinematical scale 
lengths that is currently put equal to the density scale 
length. This is partly based on Lewis and Freeman 
(1989) kinematical observations. Their kinematical gradi- 
ent seems to be conflrmed l)y few existing measurements 
including those presented in this paper. However the value 
they adopt for the density scale length is now very ques- 
tionable and is much larger than the most recent direct 
(non-kinematical) determinations. Then we cannot assess 
equality of density and kinematical scale length for our 
own Galaxy. 

If 2 7^ 0, the vertical derivative of density is not zero. 
Since we are analysing samples of stars far outside the 
mid-plane of the Galaxy, a new term appears inside the 
bracket of Equ. (2). For an exponential disk of the form 
p(z) = exp(—{z/h)"), we obtain : 

2VcVias ^ 2{A-B) ^ _ J L I In _ ^1 

+^(i--.(jr)(i-f) 

where A and B are Oort's constants and H„. H„2 are 
density and kinematical scale lengths. This relation shows 
that for an isothermal population, the asymmetric drift 
does not change with the distance above the galactic plane 
if and only if the potential is cylindrical (A = 0). We can 
estimate for typical values (A = 0.7, n=l) that the lag of 
a population with ujj = 30 km/s will vary by only 1 km/s 
on one scale height, and for ajj = 70 km/s by 5 km/s. 

6.3.1. Thick disk 

We have informations for all dominant terms present in 
the asymmetric drift relation for the thick disk popula- 
tion. We have measured the scale length directly from 
star counts, the asymmetric drift, ay in all flelds. rr^ „ 
for GC and GAC1,2 fields, and we get ujj from indepen- 
dent proper motion surveys. So we can determine one of 
these quantities from the others, and try to determine the 
self-consistency of observed quantities. 

Since the scale height of thick disk is 760 pc and the 
distances of observed stars range between 1 to 2 kpc above 
the galactic plane where the force is nearly constant 
and the density is exponential, we get n = I and z/h ^ 
1.5 — 3 in the last term of the asymmetric drift relation. We 
apply the asymmetric drift relation in order to obtain an 
estimate of ajj by replacing all the other known quantities. 
Velocity gradients are null for the thick disk population. 
Substituting V Lag = -53 km/s (with respect to the Sun), 
Vc=220 km/s, hij=2.8 kpc, R=8.09 kpc and A = 0.7 in 
Equ. (2), we obtain ajj= 72-82 km/s, in agreement with 
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Fig. 9. Asymmetric drift <V> versus (Tu,w for F and G-type stars in 2 fields (GAC1,2 and GC). The regression line gives a 
tangential motion for the Sun of Vq = 6.6±1.6 km/s 



the value given in table 5 (67 ± 4 km/s). Conversely, we 
may determine the "density" scale length by replacing ajj 
by the value given in table 5. We find from Equ. (2), the 
scale length of thick disk to be 2-2.8 kpc. 

Substitutions of all observed quantities in the asym- 
metric drift equation show the consistency of our kine- 
matical data with scale length determined by star counts. 



6.3.2. Thin disk 

In the same way, we have tried to determine the scale 
length of the thin disks from the asymmetric drift rela- 
tion. The thin disk is certainly made of a set of contin- 
uous disks with varying scale heights, lengths and kine- 
matics. This could be the reason why we do not succeed 
to extract from Equ. (2) a simple relation to predict the 
relation found on figure 9. For that purpose it would be 
necessary to know exactly the proportion of each isother- 
mal component at each height. It appears also that the 
exact value of the index n giving the shape of the vertical 
potential is determinant. Moreover, it does not seem that 
the relation should be necessary linear. Such uncertainties 
lead to very large range of possible scale length (positive 
or negative!). 



7. Discussion 

In the following subsection, we discuss various determina- 
tions of the disk's radial scale length. In the next subsec- 
tion, we discuss some of our results in terms of the history 
of the galactic disk and thick disk populations, its present 
state and the formation scenarios of the thick disk. 



7.1. The thin disk's radial scale length 

There have been several determinations of the radial scale 
length in the thin disk of the Milky Way either from pho- 
tometric or kinematic approach. 

The many different determinations based on photom- 
etry or star counts give possible values in a large range 
(Kent et al. 1991 give a recent review, see also Robin et 
al. 1992). To add to this apparent confusion, many pub- 
lished density scale lengths are deduced from kinematic 
data using the asymmetric drift relation, but generally 
the expressions used for the asymmetric drift are simpli- 
fied without any strong theoretical support. For example, 
it is frequently assumed that the kinematic and density 
scale lengths are equal. Only in a recent contribution (Fux 
& Martinet 1994) the term including the shape (spherical, 
cylindrical...) of the potential is adjusted. 

Here, we have presented new direct measurements of 
density and kinematic gradients. The resulting density 
scale length of 2.3±0.6 kpc determined by our data is the 
lower limit of currently accepted values. It has to be com- 
pared to some recent works : Robin et al. 1992, h^j = 
2.5 kpc, from star counts towards the galactic anticentre, 
Fux & Martinet 1994, Yir = 2.5 kpc, based on a rigorous 
analysis of the asymmetric drift relation. 

The kinematic gradient presented here — g^' = - 
0.18±0.03 kpc~^ is consistent with a few other direct de- 
terminations. Density and kinematic scale lengths we ob- 
tain are consistent in respect to the asymmetric drift re- 
lation. 

7.2. The form.ai'ton of the thick d,isk 

Possible thick disk formation scenarios have been sum- 
marized by Gilmore et al. (1989) and Majewski (1993). 
One of these scenarios is that the formation of the thick 
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disk occurs due to kinematically heating of thin disk. This 
heating may have been violent, perhaps from the accre- 
tion of a satellite galaxy (Carney et al. 1989, Quinn et 
al. 1993). So the merging of one or more smaller galaxies 
with our Galaxy is a possible explanation for the forma- 
tion of a thick disk. In such a scenario, the thick disk could 
be either the direct relic of the merged galaxy, dissolved 
and smoothly distributed in our Galaxy, or an old thin 
disk formed before the merger event(s) and 'puffed up' 
by the gravitational perturbations of the merging galaxy. 
In case of mergers between disks and satellites, both the 
radial and vertical structures of the disk are altered. The 
radial heating and re-arrangement of material by angu- 
lar momentum transport results in an asymmetric drift 
(Wyse 1994). The vertical heating increases the disk scale 
height by about a factor of two. Direct evidence for an on- 
going merger event has been recently discovered towards 
the galactic bulge by Ibata et al. (1994). 

Recently Quinn et al. (1993) (hereafter QHF) have pre- 
sented a model of thick disk, which is based on a scenario 
of satellite accretion by a spiral galaxy which may pro- 
duce a thick disk if the event occurs at the beginning of 
the life of the thin disk. According to this model the satel- 
lite accretions do gradually thicken disks by depositing 
ordered kinetic energy into random motions of disk stars. 
The abruptness of accretion origin implies a kinematical, 
chemical and age disjointedness from the thin disk. 

Of course the N-body simulations of QHF is only illus- 
trative, but our measurements reproduce most of features 
of their numerical simulations. Asymmetric drift (~ 50 
km/s) and scale height (760 pc) of the thick disk are in 
good agreement with the QHF model. Our measurement 
of a vertical velocity dispersion of 40 km/s for the thick 
disk is also well in agreement with the QHF accretion 
model. Also the velocity dispersion gradient of the QHF's 
thick disk is much smaller than for the disk, and is nearly 
zero in the outer part. This corresponds to the null gradi- 
ent observed. Scale lengths of thin disk and thick disk are 
also in agreement if we are in the outer part of the QHF 
model. This estimate is suggestive, however, that the thick 
disk model, which has a vertical velocity dispersion of 
45 km/s, could be formed from a thin disk with vertical 
dispersion of ~ 20 km /s by accretion of a satellite of about 
10% of the mass of the disk. 

Our measurement of the scale height of thick disk (760 
pc) is also in perfect agreement with QHF model. How- 
ever, our data do not allow to observe the correlation of 
scale height with the galacto centric radius as QHF have 
predicted in their N-body simulations. 

Finally the local surface density of the thick disk is 
about 18% of thin disk, this gives a crude estimate of the 
total mass ratio of thick disk and thin disk stars. This im- 
plies that the merging occurs fjuite early in the Galaxy life 
when 5-10% of the disk was formed. If the star formation 
rate has been nearly constant (Haywood et al. 1995) at the 
epoch of formation, this indicate that the merge occurs 2 



to 4 billions years after the formation of disk component 
(of course this primordial-disk has been transformed in 
thick disk). 

The main feature we determine for the thick disk is its 
isotliermality and its clear separation from thin disk and 
halo. This component is clearly simple and not a mixture 
of sub-components as for the thin disk and there is no 
gradient in the kinematic neither vertically nor radially. 
Similar feature has been obtained by Gilmore et al. (1995) 
and Robin et al. (1995), who find a null vertical gradient 
in metallicity : this last result can be understood as a 
consequence of the isothermality of the thick disk. 

Concerning alternative explanations, we mention 
Burkert et al. (1992) who build a 1-D (vertical) model 
of collapse of our Galaxy and formation of galactic disk. 
Their self-regulated chemical and dynamical evolution of 
an initially hot. gaseous protodisk leads to the formation 
of a thick disk. This is the only "top-down" scenario of 
formation for the thick disk where some discontinuity be- 
tween thin and thick disk is predicted with an acceptable 
agreement with our observations. However a more accu- 
rate chemico-dynamical models (Samland, 1994), with a 
2-D resolution show a more continuous age, metallicity 
and kinematic gradient between the various disk compo- 
nents. 

To sum up. the results emerging from the present study 
of the correlations between photometry and kinematics 
give a mounting evidence in favour of merging process 
of satellite galaxies with the disk of our Galaxy for the 
formation of the thick disk. 

8. Conclusion 

In this paper we have used star counts and kinematical 
data to constrain the galactic structure parameters. We 
have also shown that a proper statistical analysis of the 
data in the 5-dimensional space (V. B-V. U-B. /q. /^j) in 
comparison with the model of population synthesis allows 
one to constrain the physical model parameters. 

We obtain a value of -0.18±0.03 kpc^' for the galac- 

tic radial gradient of velocity dispersion ( — for the 
thin disk population. Determined from moderately distant 
stars, a new measurement of the solar motion is obtained. 
We found no systematic motion of the LSR relative to dis- 
tant stars on the galactic meridional plane. The rotational 
velocity curve is found fiat in the solar neighborhood. Our 
results confirm that the thin disk has a relatively short 
scale length of 2.3±0.6 kpc and scale height of 260±50 pc. 
The thick disk population is distinct from other popula- 
tions based on their kinematical and spatial distributions. 
The data constrain the asymmetric drift of the interme- 
diate population, which is found to be 53±10 km/s with 
respect to the Sun. No radial or vertical gradient is found 
in the rotational velocity and velocity dispersions of the 
thick disk population. We therefore conclude that the for- 
mation of the thick disk did not occur as a smoothly con- 
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tinuous transitional phase between formation of the halo 
and formation of the thin disk during the collapse of the 
Galaxy (Eggen, Lynden-Bell & Sandage 1962). Rather, it 
support the model of a scenario of satellite accretion by 
our Galaxy which may produce a thick disk if the event 
occurs at the beginning of the life of the thin disk (Quinn 
et al. 1993). The most probable value of scale height for 
the thick disk component is determined to be h^~760±50 
pc with a local density of A^thick = 7.4lj g % relative to 
the thin disk. The ratio of the number of thick disk stars 
in our galactic centre region to that in anticentre region 
yields hij-^Sil kpc for the scale length of thick disk. These 
values are in perfect agreement with the recent determina- 
tion given by Robin et al. (1995) based on the analysis of a 
large set of available photometric catalogues with accurate 
photometry. 
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